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^— ^ ' High-pT pion and kaon production is studied in relativistic proton-proton, proton-nucleus, and 

nucleus-nucleus collisions in a wide energy range. Cross sections are calculated based on perturbative 
^s^j ■ QCD, augmented by a phenomenological transverse momentum distribution of partons ("intrinsic 

fc-r"). An energy dependent width of the transverse momentum distribution is extracted from pion 
and charged hadron production data in proton- proton/proton-antiproton collisions. Effects of multi- 
' scattering and shadowing in the strongly interacting medium are taken into account. Enhancement 

of the transverse momentum width is introduced and parameterized to explain the Cronin effect. 
In collisions between heavy nuclei, the model over-predicts central pion production cross sections 
(more significantly at higher energies), hinting at the presence of jet quenching. Predictions are 
made for proton-nucleus and nucleus-nucleus collisions at RHIC energies. 
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PACS numbers: 24.85. +p, 13.85.Ni, 13.85.Qk, 25.75.Dw 



I. INTRODUCTION 



As the bombarding energy is increased, high-pT particle production becomes a prominent feature of nuclear col- 
1 lisions, as evidenced recently by the data beginning to emerge from the Relativistic Heavy Ion Collider (RHIC) 0. 
Qh These data call for concentrated theoretical efforts to overcome the challenges presented by the high-energy and many- 
body features of nuclear collisions at center-of-mass energies of 100 AGeV < ^/s < 200 AGeV. Relying on asymptotic 
| freedom, perturbative quantum chromodynamics (pQCD) can be applied at sufficiently high energies. To use pQCD 
, as a practical tool, one takes advantage of the factorization theorem, which provides a simple connection between the 
level of observed particles and that of the underlying quark-gluon (parton) structure. Briefly, the observable cross 
. , sections are expressed in terms of a convolution of partonic cross sections with parton distribution functions (PDFs) 
^ ■ and fragmentation functions (FFs), which encode some of the perturbatively non-calculable low-energy aspects of the 
physics. The resulting calculational scheme adopted here is referred to as the pQCD-improved parton model H[|. 
To assure the validity of a perturbative treatment, we limit the discussion to hard particle production, requiring 
the transverse momentum of the inclusively measured produced particle to be above a certain threshold, px > Pto- 
Typical values of pro ar e around 1-3 GeV in the literature. 

Hard pion and kaon production in proton-proton (pp) and proton-nucleus (pA) collisions at these energies is in itself 
of interest, but it is also imperative to study these reactions as a step in the process of understanding hard n and K 
production in nucleus- nucleus (AA) collisions. The necessity of this systematic approach has been made particularly 
clear by the high- visibility parallel developments in the subficld of J/tp production j|,U. Significant amounts of 
experimental data on ir and K production in pp collisions are available in the energy range 20 GeV < y/s < 60 GeV 
|jo|-|r2"|| . Higher energies are less thoroughly explored, with mostly calorimetric data on total hadron production (in 
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lieu of the identified pion and kaon spectra at lower energies). In the present paper, we use CERN UA1 Jl3|-p^| and 
Tevatron CDF |l6) data from pp — > /i ± reactions to bracket the RHIC energy range from above. 

The presence of nuclear effects modifies the predictions of the pQCD-improved parton model for hard particle 
production in pA collisions, as indicated by available pA hadron production data |^J^, |TT| , ^7|Jl8| . While hard particle 
production (usually disregarding nuclear effects) has been used as a testing ground for pQCD, the nuclear modifications 
are crucially important for RHIC experiments and planned nuclear experiments at the Large Hadron Collider (LHC). 
In particular, to assess the effects of jet quenching jl^-|2lj], precise background calculations and the knowledge of the 
gluon density of the medium responsible for the energy loss are necessary . 

In a recent Rapid Communication p3[ | , we suggested a specific form ( "saturation" ) of the nuclear modification for 
the cross section enhancement observed at high pt in pA collisions over what would be expected based on a simple 
scaling of the appropriate pp cross section (Cronin effect) |^|. The present paper provides a more detailed account 
of our work, incorporating not only neutral but also charged n and K production, with applications and predictions 
at RHIC energies and for AA collisions. The saturation prescription |2^] is looked upon as a limiting case, with 
the opposite limit correspondi ng to the participation of all available nucleons in the enhancement of the transverse- 
momentum width (see Section [IB ). We study the dependence of the results on the possible choices between these 
limits, and suggest a connection of the preferred value of the number of semi-hard (momentum transfer ~ 0.5 GeV) 
collisions to the results of lower-energy experiments 1 24 - 2^1. F urther measurements are needed to clarify the physical 
picture, and we argue that the planned pA runs at RHIC p8fl are likely to play a very important role in this regard. 

A similar effort (with somewhat different emphasis) was published recently [|29|| . Initial RHIC data and the interest 
in jet quenching warrant an independent detailed study of the pQCD background to which jet quenching should be 
applied to achieve agreement with central AA collision data. Our work differs from Ref p9| on the basic pQCD level 
in the choice of scales and in the use of the newly-published KKP fragmentation functions [Q. Furthermore, we 
take the position of consistently using leading order (LO) pQCD in the present paper, without introducing a so-called 
'K factor', which we found to be energy and transverse momentum dependent in another publication plfl . In this 
way we have fewer parameters, and the burden of obtaining a good fit to the pp data rests entirely on the transverse 
momentum distribution, without mixing the effects of the K factor with those of the intrinsic transverse momentum. 
We determine the energy-dependent width of the transverse momentum distribution by fitting pp data, in contrast to 
the prescribed scale dependence used in Ref. jgg ]. 

Our strategy in this paper is to first examine hard pion production in pp collisions up to yfs < 60 GeV cm. energy. 
We find that satisfactory agreement with the data can be achieved in this energy region in LO pQCD, utilizing the 
width of the intrinsic transverse-momentum distribution of partons in the colliding nucleons. We treat this quantity 
(also measured in dijet events, see e.g. Ref. |32]| ), as an energy-dependent nonperturbative parameter. This is discussed 
in Section II A , Having extracted the best value of the transverse-momentum width from the pp data, we apply these 



ideas to pion production in pA collisions in the same energy range. Our choice for the effective number of semi-hard 
collisions is described in Section II B. We then discuss AA collisions, and compare our results to CERN data on ir° 
production in S + S and Pb + Pb collisions at yfs ~ 20 AGeV. In Section [II we move to production in pp, pA, 
and AA collisions at the same ener gies, w ith the parameters as fixed above. A major step towards calculations at 



RHIC energies is covered in Section [V A , where we attempt to extend the energy range of the parameterization in 
terms of the width of the transverse-momentum distribution of partons in the nucleon. This is made difficult by the 
availability of a smaller number of pp experiments, most of them at a significantly higher y/s than the RHIC energy 
domain, limited to the measurement of total charged hadron production. The implied unc ertain ties need to be kept in 
mind when we display predictions for pA and AA collisions at RHIC energies, in Sections [VB and IV C, respectively. 
Section ^ contains our summary and conclusions. In the Appendix, we discuss the valence and sea contributions to 
quark fragmentation into charged pions and kaons. We use h = c = 1 units. 



II. PION PRODUCTION AT CENTER-OF-MASS ENERGIES BELOW ~ 60 GEV 



In this Section, we first summarize the treatment of particle production in the pQCD-improved parton model, 
listing various assumptions and ingredients. We introduce the intrinsic transverse-momentum distribution of partons 
and fit the width of this distribution to available pp data in the energy range 20 GeV < y/s < 60 GeV in Subsection 



[I A . A large body of pion and kaon production data can be utilize d in this energy interval to extract information on 
the width of the transverse momentum distribution. In Subsection [I B the Cronin effect is discussed in pA collisions. 



Subsection II C deals with hard pion production in AA collisions. 
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A. Parton model and pQCD with intrinsic transverse momentum 



The invariant cross section for the production of hadron h in a pp collision is described in the pQCD-improved 
parton model on the basis of the factorization theorem as a convolution H : 

ffii-Sr - = ^ / dx a dx b dz c fa/ P (x a ,Q 2 ) h/p(xb,Q 2 ) ^r(ab^cd) h/c ^ - sS(s + i+u) , (1) 

dp At nz£ 

c abed 



where f a / P {x, Q 2 ) and f b / P (x, Q 2 ) are the PDFs for the colliding partons a and b in the interacting protons as functions 
of momentum fraction x, at scale Q, da/dt is the hard scattering cross section of the partonic subprocess ab — > cd, 
and the FF, D h / C (z c , Q' 2 ) gives the probability for parton c to fragment into hadron h with momentum fraction z c at 
scale Q'. We use the convention that the parton- level Mandelstam variables are written with a 'hat' (like t above). 
The scales are fixed in the present work as Q = pt/2 and Q 1 — px / '(2z c ). 

In this investigation we use leading order (LO) partonic cross sections, together with LO PDFs (GRV) |3^] and 
FFs (KKP) |3(J. This ensures the consistency of the calculation. An advantage of the GRV parameterization is that 
this fit uses data down to a rather small scale (ss 0.25 GeV 2 ), and thus provides PDFs with approximate validity for 
hard processes down to relatively small transverse momenta, pt > pro = 1 GeV (with our scale fixing). It can be 
argued that, to the extent that LO and next-to-leading order (NLO) PDFs and FFs are fitted to the same data, they 
represent different parametcrizations of the same nonperturbative information. However, going to NLO will reduce 
the scale dependence of pQCD calculations pijj. A NLO study along the lines of the present work is in progress 



J55|. An alternative approach is to use the "Principle of Minimal Sensitivity" |36| to optimize these scales, as e.g. in 
Ref. [pi[ . However, fixing the scales is more convenient as a point of departure for pA and AA studies. 

Since the fragmentation functions Dh/c{ z c->Q' 2 ) are normally given for neutral hadrons or equivalently for the 
combination — (h + + hT)/2 (using isospin symmetry), assumptions need to be made to obtain the charge- 
separated FFs for 7r + , 7T~, K + , and K~ p7| . These assumptions are connected to the physical picture used to consider 
fragmentation from valence and sea quarks, respectively. The simplest approximation corresponds to allowing a meson 
to fragment only from quarks that appear in it as valence contributions. We find this prescription too restrictive, and 
use a 50-50 % division of sea quark contributions between positive and negative mesons in the present work [a = 0.5), 
unless indicated otherwise. The value of a hardly influences pion production, but has a more visible effect on kaons 
(in particular on K~). We illustrate this by also showing a = 1 results for comparison in figures with K/t: ratios. 
The details of the approximation are discussed in the Appendix. 

The PDFs in eq. (Q) express the probability of finding a parton in the proton with longitudinal momentum 
fraction x, integrated over transverse momentum up to for = Q. Recent studies have considered parton distributions 
unintegrated over the transverse momentum kr of the parton p^ j. In a more phenomenological approach, eq. (|l|) 
can be generalized to incorporate intrinsic transverse momentum by using a product assumption and extending each 
integral over the parton distribution functions to fc^-space [p9|| 



f a / p (x, Q 2 ) - > dx d 2 k T g(k T ) f a / p (x, Q 2 ) , (2) 

where g{kr) is the intrinsic transverse momentum distribution of the relevant parton in the proton. We follow the 
phenomenological approach in the present work, taking g{kr) to be a Gaussian: 

Here (k 2 -,) is the 2-dimensional width of the kj- distribution and it is related to the magnitude of the average transverse 
momentum of one parton as (k 2 ^) = A(kT) 2 /tt. In order to regularize the singularity in the cross sections (associated 
with one of the Mandelstam variables approaching zero) we use the standard procedure of introducing a regulator 
mass M — 0.8 GeV. As discussed e.g. in [ |39|| , the results in the addressed energy and transverse- momentum range 
are not sensitive to reasonable variations in this quantity. Kinematical details can be found in Ref.s |29|j39| , ^0| . 

The need for intrinsic transverse momentum in pp collisions was investigated as soon as pQCD calculations were 
applied to reproduce high-pr hadron production |2],|4l||. An average intrinsic transverse momentum of (kr) ~ 0.3 — 
0.4 GeV could be easily understood in terms of the Heisenberg uncertainty relation for partons inside the proton. 
However, a larger average transverse momentum of (fey) ~ 1 GeV was extracted from jet-jet angular distributions 
(see e.g. ]32||). Recently, new theoretical efforts were mounted to understand the physical origin of ( for) j3^,^2] , p3[ . 
Intrinsic transverse-momentum distributions have been utilized in the contexts of photon production |T8|j23| , |44| | and 



J/ip production 1 45 46J] . The description of deep inelastic scattering in the pomeron framework also appears to require 
a broad intrinsic for distribution |47| . 
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To test the validity of the above approach, we calculated ir°, ir + and ir~ production in pp collisions in the energy 
range 20 GeV < y/s < 60 GeV, and compared the results to data from several independent experiments [fr|-[l2f . The 
calculations were performed using the finite rapidity windows of the data. The Monte-Carlo integrals were carried 
out by the upgraded VEGAS-routine As discussed above, the scales are fixed as Q — pr/2 and Q' — pr/(2z c ). 
At each measured transverse momentum value (j>t > 2 GeV) we determined the width of the transverse-momentum 
distribution, (fc^), necessary to fit the data point at the given px, together with its error bar. The values of the 
extracted (fc^) are shown in Fig. 1 for 7r°, and in Fig. 2. for 7r + (top panels) and ir~ (bottom panels) production. 
The error bars are determined from the experimental errors. 
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FIG. 1. The width of the transverse momentum distribution, (fcy), necessary in eq. Q to fit the measured value of the 
spectrum at the given pr, point by point, for n° production in pp collisions at three energies. The data are from Ref. tel. 
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FIG. 2. The width of the transverse momentum distribution, (fcy), necessary in eq. (|3|) to fit the measured value of the 
spectrum at the given value of pr, point by point, for tt + (top panels) and n~ (bottom panels) production in pp collisions at 
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It is clear from Fig.-s 1 and 2 that the description of the data in terms of a pQCD calculation augmented by a 
Gaussian transverse momentum distribution makes sense only in a limited px window in the energy range 20 GeV 
< y/s < 60 GeV. In addition to the requirement of hard particle production (j>t > Pro) imposed for the applicability of 
pQCD, in most data sets there appears to be an upper limit in transverse momentum, to which (k T ) can be extracted 
with reasonable accuracy. In broad terms, the procedure may be considered sensible for, say 2 GeV < pr < 7 GeV, 
depending on the details of the experiments. Within this window, one may reasonably extract an approximately 
constant (energy-dependent) (k T ) for 3 GeV < pt < 6 GeV in the full energy range. (Note that the lower limit of the 
considered px interval can be lowered with increasing cm. energy.) At higher transverse momenta (and fixed y/s) the 
effect of intrinsic kr becomes less important, and thus it is not surprising that, using the error bars on the data, the 
width (fc|i) acquires a large uncertainty. Although 3 GeV < px ;$ 6 GeV is a rather narrow transverse-momentum 
interval, the Cronin effect in pA collisions, to be addressed in the next subsection, is most prominent in this range. 
This motivated us to use an energy-dependent, p-r-independent (fc§>) in what follows, as an approximation. 

Concentrating on the above interval in px, we next determined the ^-independent (k T ) which best fits the pion 
production data in pp collisions as a function of cm. energy. In a slight departure from our earlier work p3j, we 
fitted the data minimizing the standard x 2 — ^2(Data — Theory) 2 / '(Exp. uncertainty) 2 to obtain an optimal (fcf>) for 
each energy. Fig. 3 summarizes these results. Data are from a range of independent experiments [H-12[. 
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FIG. 3. The best fit values of (fey) in pp —> tyX J^-p^] reactions. Where larger error bars would overlap, the iv point has 
been shifted slightly to the right for better visibility. 

We conclude from this figure that the value of (k T ) for hard pion production shows an increasing trend as a function 
of y/s from 20 to 60 GeV (and may be leveling off towards the high-energy end of the interval). The numerical values 
are significantly larger than expected solely on the basis of the Heisenberg uncertainty principle. It needs to be kept 
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in mind that these values are obtained in the framework of a LO calculation. At NLO, smaller transverse width 
should be required for the best fit of the data (see some results in Ref. 31 1). At present, we focus on the utility of 
this parameterization to describe the pp data and view Fig. 3 as a basis for addressing pA and AA collisions at the 
same level of pQCD. 

As an example to illustrate the degree of accuracy of the description in the pp sector, Fig. 4 compares calculated 
7r + and 7r~ spectra and ir~ /n + ratios to the data [J7J at the cm. energies y/s = 19.4, 23.8, and 27.4 GeV, respectively, 
obtained with the values of (kj,), indicated in the top panels. We find that the data/theory and ir~ /ir + ratios are well 
reproduced for 2 GeV < pt < 6 GeV. Based on this and similar examples, we believe that hard pion production in pp 
collisions is reasonably under control at the present level of calculation. With the parameterization of the transverse 
momentum distribution we separated any pQCD uncertainties from the nuclear effects, which we address in the next 
subsection. 
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FIG. 4. Invariant cross section of ir + and 7v~ production from pp collisions (top row), D/T: data/theory ratios (middle 
row), and R: -k~/it + ratios (bottom row) as functions of transverse momentum at cm. energies y/s = 19.4, 23.8, and 27.4 GeV. 
Data are from M. The calculations are carried out in LO, as explained in the text. 



B. Hard pions from proton- nucleus collisions 



Interest in the nuclear dependence of hard particle production was triggered by the discovery of the Cronin ef- 
fect [||J7]] , and revived with the study of collisions of a particles at the CERN Interacting Storage Ring (ISR) [p)f-|5H . 
As discussed in the introduction, the nuclear enhancement of hard pion and kaon production cross sections is not only 
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interesting in itself, but we also need to understand these phenomena before we can move to the description of AA 
collisions at RHIC energies. 

The standard framework for addressing high-energy pA collisions is provided by the Glauber model [ j52|]53| ] . It should 
be kept in mind, however, that the Glauber model was originally developed with nucleons as elementary constituents, 
neglecting coherent scattering from several nucleons. Therefore, it is not surprising that refinements to the standard 
Glauber picture are proposed based on the underlying quark-gluon structure 54-56)]. As such modifications are not 
yet generally accepted or organized into a consistent picture, we continue using the standard Glauber description in 
the present work. 

The key observation for the explanation of the Cronin effect is that in pA collisions, in addition to the hard parton 
scattering, the incoming and outgoing partons may suffer additional interactions in the presence of the nuclear medium 
p9| , p9|j40| . An effective way to summarize these additional interactions is in terms of an enhancement of the width of 
the transverse momentum distribution above the value in pp collisions, shown in Fig. 3. The extra contribution to 
the width due to the nuclear environment can be related to the number of nucleon-nucleon collisions in the medium. 
To characterize the (kj,) enhancement, we write the width of the transverse momentum distribution of the partons 
in the incoming proton as 



(4) 



Here (k^) pp is the width of the transverse momentum distribution of partons in pp collisions from Subsection [I A 
(also denoted simply by (k^)), h p A(b) describes the number of effective nucleon-nucleon (NN) collisions at impact 
parameter b which impart an average transverse momentum squared C. In pA reactions, where one of the partons 
participating in the hard collision originates in a nucleon with additional NN collisions, we will use the pp width from 
Fig. 3 for one of the colliding partons and the enhanced width (Q) for the other. The function h p A(b) will be discussed 
shortly. 

According to the Glauber picture, the hard pion production cross section from pA reactions can be written as an 
integral over impact parameter b: 



d 3 P 



d 2 b t A (b) E„ 



da™((k 2 ) pA , (P T ) pp ) 



d 3 r 



(5) 



where the proton-proton cross section on the right hand side represents the cross section from eq. (|l|) with the 
transverse extension as given by eq.-s (Q) and (|J), but with the widths of these distributions as indicated. Here tAip) = 
j dz p(b, z) is the nuclear thickness function (in terms of the density distribution p) normalized as J d 2 btA(b) = A. 
Furthermore, it is well-known that the PDFs are modified in the nuclear environment ("shadowing") p?j , p8[ . We 
approximately include shadowing and the isospin asymmetry of heavy nuclei into the nuclear PDFs by considering 
the average nuclear dependence and using a scale independent parameterization of the shadowing function S a /A{x) 
adopted from Ref. 57 1: 



fa/A(x,< 



S a /A(x) 



A 



fa/p(x,Q 2 ) 



fa/n{x : 



(G) 



where f a /n(x, Q 2 ) is the PDF for the neutron. 

The effectivity function h p A(b) can be written in terms of the number of collisions suffered by the incoming proton 
in the target nucleus, v A ip) = ONNtA{b), where <Jnn is the inelastic nucleon-nucleon cross section. Two extreme 
prescriptions were used for h p A(b) earlier: (i) h^^b) — VA{b) — 1, corresponding to all NN collisions, except the hard 
interaction producing the parton to fragment into the observed pion, contributing to the (fc^) enhancement fl29|] , and 
(ii) a 'saturated' prescription, where hp^fb) was equated to a maximum value of unity whenever VA{b) > v m = 2 (i.e. 
it was assumed that only one associated NN collision is responsible for the (k^) enhancement, further collisions are 
not important in this regard) p3| . In the present work, we examine the dependence of the results on the possible 
choices between these limits, denoted by v m = oo and v m = 2, respectively. Since even in a heavy nucleus VA{b) does 
not exceed « 6, we carried out explicit calculations for 2 < v m < 5 and without any cut (v m — oo). The coefficient C 
in eq. ([|) best describing the available pA data for the different values of v m is shown at three energies in Fig. 5 
as a function of px for Be, Ti, and W target nuclei. 
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FIG. 5. The best fit values of C in eq. (^) for Be, Ti, and W targets at three different energies as a function of pr and u„ 
The fitted data are from Ref. [q|. 



For an appealing physical interpretation, the coefficient C is expected to be approximately independent of pt, 
of the target used, and probably of beam energy (at least in the energy range studied in Fig. ||). Based on these 
requirements, we select v m — 4 as the value that produces results for C closest to the ideal physical picture. In what 
follows, we fix v m = 4, where the associated value of C is C = 0.40 ± 0.05 GeV 2 . This value of C is significantly 
smaller than what we needed to use in connection to v m = 2 earlier (C sat , see Ref. [^3|), and somewhat larger than 
would be needed in the case of v m = oo. This can be understood based on the compensating effects of h p A{b) and 
the value of C in eq. (||). 

It is interesting to note that a similar value for the number of effective inelastic (" semi- hard" ) collisions was obtained 
examining nuclear stopping in p + Be,Cu, Au collisions at lower energies (E}, eam =12 and 18 AGeV) p^-|2^|. The 
leading proton appeared in these experiments not to loose more energy after 2-3 collisions (measured via 'grey' tracks 
from the target). The physics of this phenomenon was connected to the constituent quark picture of the proton, i.e. it 
was suggested that the state resulting after the inelastic excitation of the three constituent quarks no longer interacts 
with a high probability. Further experiments are needed to verify this or a similar physical interpretation. Such an 
effort is in progress at AGS and at CERN SPS (see Ref. [^6|,^7)). Our finding of v rn — 4, which means ~ 3 effective 
semi-hard collisions, overlaps with the above result and points to a similar mechanism, even at higher energies. At 
this point, the present study lends phenomenological support to the explanation advanced in Ref. p4|-p6|]. Our model 
is also consistent with the idea proposed by Dokshitzer |59] about the proton being unable to accept more than a 
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characteristic momentum transfer of Q w 1 — 1.2 GeV. This idea implies the existence of a limiting value of (fcf,) for 
the proton, which can be accumulated in one or several semi-hard interactions. Further data are needed at CERN 
SPS and RHIC energies for a more definitive statement. The planned pA collision program at RHIC p8| may play a 
very important role in this regard. 

To illustrate the quality of agreement with the data achieved with the above parameterization, Fig. 6 displays 
7r + spectra (upper panels) and ir~ /n + ratios (lower panels) as functions of px, compared to the same pA data jjj 
as used to determine C. The parameters are fixed as indicated in the top panels: v m = 4, C — 0.4 GeV 2 , and the 
energy-dependent (k T ) pp from Fig. ||. There is no observable difference in the ir~ /tt + ratios between the different 
targets in either the data or the calculated results. 
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The Cronin enhancement can be presented in the form of normalized (by mass number A) W/ Be cross section ratios 



in Fig 



The data for tt + and 7r~ are shown as a function of px for the above energies and for Ei a b — 800 GeV |11| ] 
together with the results of calculations with C = 0.35 GeV 2 (solid line) and with C — 0.45 GeV 2 (dashed 
line). The band between these two curves represents the uncertainty in our calculation associated with the extraction 
of the coefficient C (see Fig. ||). 

In the absence of the Cronin effect, these ratios would give a value of identically 1. The significant deviation of the 
data from unity is therefore a clear confirmation of the nuclear enhancement in the 2 GeV < pr < 6 GeV transverse 
momentum window. At lower pt, where absorption effects set in, the data indicate a value of the ratio smaller than 
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one. However, we do not trust our pQCD calculation below px ~ 2 GeV, and thus do not show calculated results 
for low transverse momenta. At high transverse momenta, where the effects of intrinsic kx become unimportant as 
mentioned earlier, the ratio converges to unity in both the data and the calculations. 
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FIG. 7. Normalized W/Be charged pion cross section ratios at four energies. The calculations are carried out with the 
parameters fixed earlier as discussed in the text. The data are from Refs. J^, |lo|Jlll . The deviation from unity represents the 
Cronin enhancement. 

C. Nucleus- nucleus collisions 

In nucleus- nucleus reactions^, where both partons entering the hard collision originate in nucleons with additional 
semi-hard collisions, we apply the enhanced width of the parton distribution (01) for both initial partons. Thus, 



E„- 



da 



AB 



d 3 ? 



d 2 bd 2 r t A (r) t B (\b-f\)E„ 



d 3 p 



(7) 



where the proton-proton cross section on the right hand side represents the cross section from eq. (|l|) with the 
transverse extension as given by eq.-s (g) and (||), with the enhanced widths of these distributions given by (|J). We 
emphasize that the parameter values used to describe hard pion production in AA reactions are identical to the ones 
used earlier in this Section for the best description of pion production from pp and pA reactions. In other words, the 
values of (fcy) pp =1.6 and 1.7 GeV 2 , respectively, are taken from Fig. || (with a relatively small uncertainty of ±0.1 
GeV at y/s ~ 17- 20 GeV), and v m = 4 and C = 0.4 GeV 2 are fixed. 

The results of these calculations for central tt° spectra are compared to the WA80 |6^| and WA98 |3l| data for 
S + S, S + Au, and Pb + Pb collisions in Fig. Has a function of transverse momentum. Spectra are shown in the top 



1 These are denoted by AA in the text, even though not all studied systems are symmetric. 
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panel, data/theory ratios are displayed at the bottom. The dotted lines represent the results of the pQCD calculation 
with nuclear effects (shadowing and multiscattering) turned off; solid lines correspond to the full calculation. When 
comparing to the data, it should be kept in mind that the calculation is intended for hard particle production, and 
it is not expected to describe the data below pt ~ 2 GeV. At low pr, where soft processes become important, the 
pQCD model underestimates the S + S data, as expected. The calculation without nuclear effects (dotted lines) 
underestimates the data in the entire measured px range. The reproduction of the S + S and S + Au data can be 
considered reasonable for pr > 2 GeV, while the Pb + Pb calculation overestimates the data in the same transverse 
momentum window by upto 40%. This may be taken as a hint that an additional mechanism is at work in the nuclear 
medium, which acts to reduce the calculated cross sections. Jet quenching p~0|— [2l]] is a potential candidate for the 
physics not included in the present model. 




Pt (GeV) 

FIG. 8. Neutral pion production compared to data from SPS experiments WA80 [[so) and WA98 for central collisions. 
Calculated invariant cross sections with the nuclear effects turned off (dotted lines) and with the complete model (full lines) are 
displayed on top together with the corresponding data, data/theory ratios are shown on the bottom (with the full calculation). 
The parameters are fixed by pp and pA information as discussed in the text. 



III. HARD KAON PRODUCTION AT CM. ENERGIES BELOW ~ 60 GEV 



In this Section we turn to charged kaon production in the same cm. energy region, using the parameter values 
fixed in Section |l[ We are interested in the reproduction of kaon data by including the effects of multiscattering and 
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shadowing. Most data are available as kaon to pion ratios; we will display the calculated results in the same manner. 
Furthermore, as mentioned above, K~ production is particularly sensitive to assumptions made in connection with 
sea quark fragmentation, discussed in the Appendix. Our standard computations are carried out with the value of 
(7 = 0.5 for the parameter governing the distribution of sea quarks, but as a reminder to this dependence and to 
illustrate the sensitivity, we also show calculated results with a = 1 in this Section. 



A. Charged kaons from proton-proton collisions 

First we look at kaon spectra and K/ir ratios for positive and negative mesons from pp collisions to check if the 
parameterization introduced for pions is also applicable in this case. Our results are compared to the data in Fig. ^ 
for K + and in Fig. [l(] for K~ at the cm. energies y/s — 19.4, 23.8, and 27.4 GeV. In the top panels the appropriate 
pion spectra are also shown as a reference. The kaon cross sections displayed here are calculated with the same value 
of a = 0.5. The data/theory ratios show an agreement similar to that seen in Fig. |]for pions. In the bottom panels 
presenting the K/ir ratios, full lines correspond to the standard a — 0.5, dashed lines depict the calculated results 
with cr = l. 
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FIG. 9. Spectra of tt+ and K A 
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mesons as functions of transverse momentum from pp collisions at cm. energies y/s — 19.4, 

0.5 (solid) and a = 1 (dashed) as explained in the text 
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FIG. 10. Spectra of 7r~ and K~ mesons as functions of transverse momentum from pp collisions at cm. energies yfs = 
19.4, 23.8, and 27.4 GeV (top panels) and K~ /ir~ ratios calculated with a = 0.5 (solid) and a — 1 (dashed) as explained in 
the text (bottom panels). Data are from |7J. 



The reproduction of the kaon data from pp collisions can be considered fair, except that in the K/tt ratios the 
apparent structure of the K + /it + ratios is not reproduced, and the peak in the K~ /tt~ ratios is shifted to higher 
transverse momenta and is less pronounced than in the data. It can also be seen that the distribution of sea quarks 
for the fragmentation of charged mesons, while not too important for the K + /ir + ratio, makes a significant difference 
for K~ /ir~ , with a — 0.5 describing these ratios better than a = 1 for pt < 5 GeV. If we focus on pt > 5 GeV, the 
opposite conclusion may be drawn. The sensitivity of the negative ratios to a is also larger in the E605 experimental 
K/tt ratios shown in Fig. [ll]. Here, measurements extend to higher values of pr, and a = 1 gives a better fit to 
the K~ /n~ ratios for p^ > 5 GeV, similarly to the K~ /n~ ratios of Ref. [Q. This may be taken as a hint that our 
prescription for the fragmentation contribution of sea quarks is too schematic, and a should be taken p^-dependent 
in a more realistic calculation. 
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FIG. 11. Ratios of K + /n + (top) and K~ /-k~ (bottom) at Ebeam = 800 GeV, calculated with a — 0.5 (solid) and a — 1 
(dashed) . Data are from Jl2,|ll| . 



B. Kaons from proton- nucleus collisions 



Next we examine the predictive power of our model for charged kaon production from pA collisions at cm. energies 
yfs = 19.4, 2 3.8, an d 27.4 GeV, with the parameters fixed as above. In other words, (fef.) is taken from Fig. |^ as 
in Subsection III A , and the enhancement of the transverse momentum width is given by eq. (Q) using u m = 4 and 
C = 0.4 GeV 2 . 

The data are available as K/tt ratios, eliminating experimental normalization issues |Q. In Fig. [l^ we compare 
the results of our calculations to these data for Be,Ti, and W targets. K + /ir + ratios are shown in the top panels, 
while K~ /ir~ ratios are displayed on the bottom. We present calculations with both a = 0.5 and a = 1. The general 
tendency of these results and the quality of the agreement between the data and the calculations is very similar to that 
observed in pp collisions. The deviations increase with energy, but the calculations seem to reflect the gross features 
of the data. The kaon to pion ratio is not sensitive to the target size. This means that the nuclear enhancement is 
similar to pions and kaons. 
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FIG. 12. K + /tv + (top) and K ~/tt~ (bottom) ratios for three nuclei at cm. energies y/s = 19.4, 23.8, and 27.4 GeV. The 
two bands of calculations correspond to a = 0.5 and a = 1, respectively. Data are from |7J. 



C. Predictions for kaons from nucleus- nucleus collisions 

We are not aware of hard (pr > 2 GeV) kaon data from AA collisions at SPS energies. However, it is natural 
to extend the present calculations to kaon production under the conditions of the WA80 and WA98 experiments. 
Predicted K + /n + (top) and K~ /n~ (bottom) ratios for S + S collisions at y/s = 19.4 GeV (left panels), and Pb + Pb 
collisions at y/s — 17.3 GeV (right panels), as functions of px are displayed in Fig. [l^. 

The appearance of these results is similar to that of the calculated results in the pA case, except that the curves 
are less steep for AA collisions. Any possible peak at transverse momenta below p? = 2 GeV in the K ~ /n~ ratios 
would be unaccessible for our calculation. Information on these ratios can serve as a test of the present model at SPS 
energies. 
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FIG. 13. Predicted K + /ty+ (top) and K~ /ty~ (bottom) ratios for WA80 (S + S at sfs = 19.4 GeV, left panels), and WA98 
(Pb + Pb at y/s = 17.3 GeV, right panels) experimental conditions. The parameters of the calculation are fixed as discussed in 
the text. Full lines correspond to a = 0.5, dashed lines display results with <x = 1. 



IV. EXTENSION TO RHIC ENERGIES 



In previous sections, we have developed a description of hard pion and kaon production based on the pQCD- 
improved parton model, incorporating a phenomenological treatment of the transverse momentum degree of freedom 
of the partons. We have seen that a satisfactory agreement with pp data can be achieved up to y/s as 60 GeV, 
allowing the width of the transverse-momentum distribution, (kj), to depend on energy. The available significant 
body of information on hard pion production in the energy range 20 GeV < y/s < 60 GeV strongly constrains the 
value of (fe|,) at lower energies, but the uncertainty increases with y/s. Modeling the enhancement of the transverse- 
momentum width associated with additional collisions in the nuclear medium, we were able to reasonably reproduce 
hard pion-production data from pA and AA reactions in the same energy interval. We then illustrated the use of the 
model for kaons. 
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A. Hadron production in pp collisions at energies ^J~s > 100 GeV 



In accordance with the main motivation for this study in terms of emerging data and future experiments at RHIC, 
in this Section the parameterization of the pp data will be extended to the RHIC energy region. Unfortunately, we 
are not aware of useful identified hard meson production data from pp or pp reactions at energies above y/s ~ 60 GeV; 
we have total hadron production data in the form of (h + + h~)/2 from CERN UA1 ]T^-[l5t and the Tevatron CDF 
fiiif . Since the production of hadrons other than kaons and pions (in particular protons and antiprotons) increases 
with energy, the accuracy of the extraction of (fcy) decreases with increasing energy. Furthermore, we have to face a 
shortage of data in the energy interval 60 GeV < ^/s < 500 GeV. 

In Fig. [l4| we show the quality of the fit we can achieve with the energy-dependent (fc^) to the (h + + hr)j2 data 
from pp collisions at ^fs =200, 500, and 900 GeV juj. The top panel contains the spectra, the bottom panel gives 
the data/theory ratios as a function of pt for the values of (/c|>) indicated in the top panel. 
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FIG. 14. Top panel: spectra of charged hadrons (h + + h~ )/2 from pp collisions at ^fs =200, 500, and 900 GeV as calculated 
with the indicated values of (fcy). Bottom panel: data/theory ratios. The data are from the UA1 experiment Jl4| , 

Fig. |l5| summarizes the extracted values of the transverse momentum width from the {h + + hr)j2 data of the 
higher-energy experiments |13|-|l6|, together with the points at lower energies extracted from hard pion production 
data (shown earlier in Fig. a). The solid curve is drawn to guide the eye, and is an attempt to represent the energy 
dependence of the extracted (fcy). The dashed lines delineate an estimated band of increasing uncertainty with 
increasing ^/s around this curve. We will use this band to interpolate between the low-energy and high-energy parts 
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of the Figure, reading off the (k^) pp values to be applied at RHIC. The darker bars represent an estimate at around 
y/s = 130 and 200 GeV, respectively. 
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shifted slightly for better visibility. The band is drawn to guide the eye. 

It can be seen that (kj) reaches a maximum at around y/s = 60-100 GeV, with a value of « 2.4 GeV 2 . Beyond 
this energy, the transverse-momentum width necessary to reproduce hadron production data from pp or pp reactions 
at the LO level of the pQCD calculation starts to decrease. The interpolating curve is drawn so that it converges to 
a small constant value of {kj) as yfs — > oo. The high energy points are based on charged pion and kaon production 
only, and should be looked upon as upper limits. In the RHIC energy region, we extract (fc|,) = 2.0 ± 0.4 at i/s w 
130 GeV and (kj,) = 1-7 ± 0.4 at y/s w 200 GeV for the value of the transverse-momentum width of partons without 
nuclear enhancement. 

For the pA and AA calculations we fixed v m = 4 and C = 0.4 GeV 2 , in lack of a less arbitrary prescription. It 
needs to be kept in mind, however, that these quantities and the efiectivity function h p A(b) of Eq. |] may well be 
energy dependent. We carried out computations for pion and kaon production at the RHIC energies of the 2000 run 
and of the 2001 run of ^/s =130 GeV and 200 GeV, respectively. 
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B. Hadron production in pA collisions at RHIC energies 



We are now prepared to make certain predictions concerning the planned pA program at RHIC pq ]. Since the 
precise energies and targets to be used are not pre sently available, we use the cm. energies of ^/s = 130 and 200 GeV, 
respectively, and the targets familiar from Section II B, as an illustration. For the width of the transverse momentum 
distribution of partons in the proton, we take the values (fc^) = 2.0 GeV 2 and 1.7 GeV 2 at y/s =130 and 200 GeV, 



respectively, representing the centers of the dark bars in Fig. 15. The parameters v m = 4 and C = 0.4 GeV 2 are 
fixed. In the top panels of Fig. |l6| we show calculated invariant 7r° production cross sections at these energies from 
collisions of protons with Be (solid line), Ti (dotted), and W (dot-dash) targets. We also include the results of pp 
calculations as a reference, and the UA1 pp data at ^/s = 200 GeV for production for comparison |L4]]. These data 
were used earlier in Fig. |l4j. Since kaon production typically does not amount to more than 20% of pion production 
and the use of antiprotons is not expected to make a significant difference at these energies, the p + p — > tt° + X 
calculation appears to fit these data on the logarithmic scale almost as well as the earlier /i ± calculation (see Fig. [Ti]). 
In the bottom panels of Fig. [l(] K + /t: + ratios are displayed for the p + W reactions, at the values of a — 0.5 (solid 
line) and a = 1.0 (dashed) for the parameter controlling sea-quark fragmentation. These ratios are very similar for 
the other targets. 
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FIG. 16. Top panels: predicted invariant cross sections for ty° production from p + Be/Ti, W collisions at yfs = 130 GeV 
(left) and 200 GeV (right). Calculated pp results and the UA1 h data from pp collisions [n4j are also included. Bottom panels: 
predicted K + /ir + ratios for p + W collisions with a — 0.5(solid) and a — 1.0 (dashed). 

Figure [l^ displays similar predictions for K + production from pA collisions. In the bottom panels we now show 
K~ /K + ratios in p + W collisions for the two values of the parameter a used earlier. It can be seen from the bottom 
panels that the value of a plays a more important role in kaon than in pion production also at these higher energies. 
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FIG. 17. Top panels: predicted invariant cross sections for K production from p + Be,Ti,W collisions at ^/s — 130 
GeV (left) and 200 GeV (right). Calculated pp results are also included. Bottom panels: predicted K~ /K + ratios for p + W 
collisions with a = 0.5(solid) and a = 1.0 (dashed). 



C. Hadron production in Au + Au collisions at RHIC 

In the area of AA collisions, calculated high-p^ ^ production in Au + Au collisions is compared to preliminary 
PHENIX data at yfs =130 GeV in the top left panel of Fig. |lj; predictions at 200 GeV are displayed on the right 

2.0 GeV 2 and 1.7 GeV 2 at ^/s =130 and 200 GeV, respectively, representing 

T2 



of Fig. 18, We use the values 



the centers of the dark bars in Fig. tLB. The parameters v m — 4 and C — 0.4 GeV 2 are fixed. Central and peripheral 
collisions are defined by the top 10%and the 60-80% bin of the total cross section, corresponding to the experimental 
selection. We also show the results of the pQCD calculations of hard pion production from pp collisions for reference. 
The pp data on total charged hadron production at 200 GeV are also included in the top right panel of the Figure 
(symbols). In the bottom panels the predicted K + /tt + ratios associated with central collisions are displayed. The 
corresponding peripheral ratios are almost identical to the central ratios shown. 
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FIG. 18. Top left: calculated central 10% (full line) and peripheral 60-80% (dotted) hard ty° spectra from Au + Au collisions 
at y/s =130 GeV, compared to preliminary PHENIX data We also show the result of our pp calculation (dashed). Top 
right: prediction for central (full line) and peripheral (dotted) hard n° spectra from Au + Au collisions at yfs =200 GeV. The 
pp results are included and compared to the UA1 pp data at the same energy Bottom: predicted K + /ir + ratios at the 

same energies. 

It can be seen in the top left panel of Fig. [l8| that, while our pQCD calculation augmented by nuclear effects ap- 
proximately reproduces the data on peripheral Au + Au collisions at \fs =130 GeV, the central data are overestimated 
by a factor of 3 — 5. This indicates that an additional mechanism is at work in dense nuclear matter, which would 
decrease the calculated cross sections at a given px- Since this effect can also be looked upon as a shift of the spectra 
to lower pt, we speculate that the phenomenon of jet quenching [Ts|-[2l|] is responsible for the discrepancy in central 
collisions. A study of jet quenching in the present framework is in progress pq]. 
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FIG. 19. Predicted K + spectra and for K~ / K + ratios for Au + Au reactions at yfs = 130 GeV (left panels) and 200 GeV 
(right panels). 



The predictions in Sections IV B and IV C above are intended to illustrate the capabilities of the description we 
developed. Similar results can be obtained easily for the precise conditions of future RHIC measurements. We find a 
systematic study of the energy dependence of hard pion and kaon production in pp collisions particularly important. 
Similarly, a study of the energy and target dependence of pion and kaon production in pA collisions will help pin 
down the energy dependence of the parameters that characterize the enhancement of the transverse momentum width 
in the medium, which have been taken energy independent for the time being. 



V. SUMMARY AND CONCLUSION 



The commissioning of RHIC opened an exciting new era in nuclear collision physics, with the study of excited 
strongly-interacting matter becoming a reality. It is recognized, however, that a thorough understanding of ultrarel- 
ativistic AA collisions presupposes the accurate description of pp and pA collisions in the same framework. In the 
present paper we attempted to follow the evolution of hard pion and kaon production in relativistic collisions from pp 
to pA to AA reactions. 

The pQCD-improved parton model suggested itself as a natural tool for our study. Of course, pQCD itself is 
evolving from leading-order calculations through the increasing complexity of NLO and NNLO for selected processes. 
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As our major focus in this paper was on the additional physics brought in by pA and AA collisions, it was decided 
to use LO pQCD throughout. This was supported by evidence that higher order pQCD does not eliminate the 
need for the transverse momentum distribution of partons for a satisfactory fit of pion and kaon production data in 
pp collisions in the 2 < px < 6 GeV window (though the numerical values become smaller, as expected). We used 
abundant pion production data from pp collisions at cm. energies y/s < 60 GeV to extract the width of the transverse 
momentum distribution of partons in the nucleon. The phenomenological value of the description was then tested on 
kaon production at the same energies. 

For the treatment of nuclear systems, we developed a model based on the enhancement of the width of the transverse 
momentum distribution of partons in the nuclear medium. An additional parameter was fitted to describe the Cronin 
effect at these energies. Shadowing and the isospin asymmetry of heavy nuclei were taken into account. We tested the 
model on charged pion and kaon production. In AA collisions at SPS energies we found an indication of a possible 
need for an additional mechanism to decrease the calculated cross section of neutral pion production in the collision 
of heavy nuclei. We speculated that jet quenching may provide that mechanism, beginning to appear already at SPS 
energies. 

Using pp data at higher energies, the domain of the extracted transverse momentum width was extended to cover 
planned and present RHIC energies. Keeping in mind the rather large uncertainties, predictions were made for pA 
and AA collisions at RHIC. The over-prediction of the preliminary PHENIX data on central tt° production at y/s = 
130 GeV reinforced our tentative conclusion that additional physics needs to be incorporated in the model to decrease 
the calculated cross sections when dense nuclear matter is present. Since the disagreement with the central data is 
much stronger at these energies than at SPS, and jet quenching is expected to increase with energy, jet quenching is 
a likely candidate for this physics. 

Recognizing that the treatment of jet energy loss would require further assumptions and modeling, we leave it for 
later study. A future treatment of jet quenching in the same framework promises to be of interest as a measure of 
the gluon density of the medium responsible for the energy loss. We are also interested in a similar study at the NLO 
level of pQCD. 

In conclusion, it appears that a description of pA and AA collisions based on the pQCD-improved parton model 
augmented by the transverse momentum distribution of partons and by the nuclear effects of multiscattering and 
shadowing works reasonably well as a background calculation for hard pion and kaon production at sufficiently high 
energies and has useful predictive power for the interpretation of present and future RHIC data. A full treatment 
of nuclear effects will need to take jet energy loss into account. We see our model as a flexible tool to aid in the 
understanding of the properties of extended strongly- interacting matter created at RHIC, LHC, and potential future 
nuclear colliders. 
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VII. APPENDIX: FRAGMENTATION FUNCTIONS FOR CHARGE-SEPARATED PIONS AND KAONS 

The parameterization of fragmentation functions (FFs) for pions and kaons is usually given in terms of neutral 
linear combinations of the charged mesons Uf0|, denoted by e.g. DJ ± , DJ*, and D"f for pions. Since the present 
paper focuses attention on separately measured 7r + , 7t~, K + , and K~ spectra and their ratios, we need separate FFs 
for the charged pions and kaons. In this Appendix we describe our approximation used for the charge-separated pion 
and kaon FFs. 

To lowest order, one may assume that each meson species will only fragment from those quarks that appear in it 
as valence quarks, or from gluons. In this approximation, there is e.g. no contribution to the fragmentation of pions 
from strange quarks, or to the fragmentation of kaons from d or d quarks. However, these flavors can appear as sea 
quarks in the respective mesons, and we want to take into account the possibility that pions or kaons fragment from 
quarks that appear in them as sea contributions. This implies a separation of the FFs into valence and sea parts, 



23 



along the lines of Ref. |3ffl . For consistency, the sea quarks corresponding to the valence flavors in the given meson 
need also be considered. 

Based on the above ideas, we split the fragmentation functions for pions as 

D u = D u,val + D Z,sea ( 8 ) 
Df = D tal + Dtsea (9) 
Df = D^ sea ■ (10) 

The sea contribution from the fragmentation of a it or d quark can be identified in this case with the FF of the s 
quark, which appears in the pion only as a sea quark, 

^u,sea -^d.sea -^s.sea ' (H) 

Furthermore, it is natural to make the usual ansatz connecting the charge-averaged FFs of quarks and antiquarks 
&■ 

Df = Df , (12) 

and similarly for d and s quarks and antiquarks. 

The valence contributions to the tt^ FFs can be obtained from the parameterization available in Ref. |3(| as 

DlU = - Df (13) 

(14) 

Since u and d occur in the positive pion as valence quarks, we have 

^u,val = ^u,val (15) 

m + , = m± n 6 ) 

d,val d,vaL 1 \ ' 

and similarly for 7r~ . 

The remaining task is to divide the contribution of u quark fragmentation to pions between 7r + and 7r _ when they 
appear in the created pion as sea quarks. We propose 

^u,sea ® ^u,sea (I'O 

Ksea = (1 - °)V>t S e a (18) 

Ksea = (1 - <r)Dtsea (19) 

D l,sea = a D l,sea ( 20 ) 

where < a < 1 is a free parameter. It is natural to expect a = 0.5 from symmetry arguments. We find that a = 0.5 
indeed provides a satisfactory description of the available pion data, and use this as the default value in the present 
paper. 

For charged kaons a very similar analysis, with the s quark replacing the d quark as a valence contribution and the 
d quark playing the role of the contribution that can only arise from the sea, defines the fragmentation functions in 
an analogous manner. 



u,val 


^ Li 


-Ds 


)t 

d,val 


— 7T 7r± 

— u d 
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